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Abstract 

We present a new method for calculating ground state properties of quantum dots in high 
magnetic fields. It takes into account the equilibrium positions of electrons in a Wigner 
cluster to minimize the interaction energy in the high field limit. Assuming perfect spin 
alignment the many-body trial function is a single Slater determinant of overlapping oscil- 
lator functions from the lowest Landau level centered at and near the classical equilibrium 
positions. We obtain an analytic expression for the ground state energy and present numer- 
ical results for up to = 40. 

1. Introduction 

Quantum dots, often referred to as artificial atoms []T|, were intensively studied in recent 
years both experimentally 0,^ and theoretically |^^18|. A number of different theoretical 



methods, including, e.g., exact diagonalization quantum Monte Carlo (QMC) ^ 



12| , Hartree-Fock ||13Hl5| and density-functional calculations ||T6Hl8|, were employed for 
investigating the ground state properties of dots, which are responsible for many of the 
experimentally observed effects. 

In this paper we present a new concept to calculate ground state properties of quantum 
dots in high magnetic fields, with application to systems with up to 40 electrons. This 
approach is motivated by the evidence that in high magnetic fields B quantum dot electrons 
tend to form a Wigner cluster Therefore, we construct a ground-state many-body 

variational wave function in the form of a Slater determinant with (overlapping) single- 
particle wave functions centered at or near equilibrium positions of classical point particles. 
The classical configurations are calculated either analytically (for A^ < 8) or using the Monte 
Carlo minimization |]8|JT9[| . Expectation values of observables can be formulated analytically. 



Here we present ground state energies for disc-like quantum dots with up to 40 electrons. 

2. The model Hamiltonian 

We consider a system of A^ electrons moving in the plane z = 0, in the lateral potential 
V{r) = rrfujl^x^ + ?/^)/2, and in the perpendicular magnetic field B = (0,0,5). The 
Hamiltonian of the system has the form 



i=l 



p. + eA(r.)f ^^^^^ 



2m* 



+ Vc + E2 
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where A(r) = {B/2){—y,x), Vc is the Coulomb energy of electron-electron interaction and 
Ez is the Zeeman energy. The solutions of the single-particle problem (Fock and Dar- 
win, [|0l) are given by ipn,i{r,ip) oc e'''^rl'le"''^/^L.ll'(r^), where = (x^ + y^)/\^, = 

the cyclotron frequency, L„ are associated La- 
guerre polynomials, n G No is the number of nodes, and / e Z is the angular momentum 
quantum number. 

3. The new approach 

In strong magnetic fields electrons tend to localize around the classical equilibrium posi- 
tions, and their spins to be aligned along the magnetic field. This tendency suggests to write 
the trial many-body wave function ^ i in the form of a Slater determinant (as has been done 



in Ref. [21 1 for a related problem) 



= det l/v^, xlf = ^o,-L(r. - R,) exp (-^ (B x R,) ■ r,^ , (2) 
where xif^ are Fock-Darwin wave functions with an additional phase factor caused by a 



magnetic translation |^ of the i-th electron to the center at R^. The non-negative integer 



L and the vectors Ri, . . . , R^v are free parameters of the theory. In general, the single particle 
wave functions Xi^^ are not orthogonal, and the overlap of the neighbor single-particle states 
needs not to be small. 

Using the trial wave functions we analytically calculate expectation values of the 
energy and the density of electrons. The energy is minimized with respect to the variational 
parameters L and Ri, . . . , R^v- 

4. Results 

Figure ^ shows the calculated ground state energy of a system of 3-10 electrons, together 
with reference data from Ref. [§ , obtained by the QMC method. The small numbers and 
1 over the dotted lines indicate that, starting from the corresponding points the trial wave 
functions \&l=o or \I/l=i gave the lowest energy. The parameters Rj were chosen according 
to classical equilibrium configurations and were not varied in Figure |I|. 

As seen from the Figure, our results are very close to (and sometimes even lower than) 
the QMC results at uJc/ujo > 2. In this range typical relative deviations of the total energy 
are 2% and less. Optimizing the parameters Rj we could get even better agreement with 
the QMC results in the range 0.6 < uJc/ujo < 2.4 (the optimized radius of electron shells is 
up to 20% larger at small magnetic fields than the classical radius). The accuracy of our 
results shows a tendency to increase with larger B and/or N. At the upper range of the 
magnetic field the relative deviation is about 1% for A < 6 and even less for higher A. This 
remarkably good agreement and its tendency to become even better with increasing electron 
number demonstrates the capability of the new method and is promise for obtaining reliable 
results for even higher electron numbers. 

Figure ^ gives calculated ground state energies in the range of up to A = 40, for which 
suitable exact reference values were not available. With increasing magnetic field we find a 
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Figure 1. Calculated expectation 
values for the ground state energy 
(dotted) of = 3 — 10 electrons. 
Bold numbers designate the num- 
ber of electrons, u = {lJq + a;^/4)^^^. 
Full curves show the QMC results 
from Ref. 0. For this plot we 
have chosen lo/as = 1.91, where 
/q = h/m*ujQ and as is the effective 
Bohr radius. The Zeeman energy 
with effective g-factor g* = —0.44 
was taken into account. 
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Figure 2. Calculated ground state 
energies for 10-40 electrons (dotted 
curves) . The horizontal dashed lines 
indicate the value of classical ground 
state energy E'^^"''^^ / hjjQ. The mean- 
ing of the numbers in the plot and 
the parameter lo/as are as in figure 
ID, except g* = here. 



clear convergence of Ej^ — Nhw towards ii^^"**, so that we can write 



E 



N 



(3) 



The correction 5Em is positive, tends to zero at i? ^ oo, and for 10 < < 40 it is smaller 
than 2% at ud^o > 6 and smaller than 0.6% at Uc/uJo > 18. Equation (H) refines the rigorous 



mathematical result by Yngvason and has a transparent physical meaning: The quantum 
mechanical ground state energy of the dot is a sum of the classical potential energy, the 
quantum mechanical kinetic energy, and the small rest energy ^-Eat, which decreases with B. 



5. Conclusion 

We have presented a new method to calculate ground state expectation values for parabolic 
quantum dots in (high) magnetic fields. The results were obtained by evaluating the expec- 
tation value, for which we found an analytical expression using a specially designed trial wave 
function. It takes into account the classical equilibrium positions of the electrons and has 
the correct symmetry for a fermion system with fully polarized spins. In the high magnetic 
field regime our results compare favourably with QMC data, but are much less demanding 
in computer time. Therefore, an extension to higher electron numbers becomes possible as 
we demonstrate by showing results for up to 40 electrons. We put forward a simple inter- 
pretation of our results for the ground state energy, Eq. (^, and gave the upper estimate 
for the rest energy term SEj^i in the studied range of N and the magnetic field. 

This work was supported by the Deutsche Forschungsgemeinschaft within the Graduier- 
tenkolleg (GK 176) Komplexitat in Festkorpern. 
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